Slow spin-glass and fast spin-liquid components in quasi-two-dimensional La2(Cu,iv^)04 
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In conventional spin glasses, magnetic interaction is not strongly anisotropic and the entire spin 
system is believed to be frozen below the spin-glass transition temperature. In La2Cuo.94Lio.06O4, for 
which the in-plane exchange interaction dominates the interplane one, only a fraction of spins with 
antiferromagnetic correlations extending to neighboring planes become spin-glass. The remaining 
spins with only in-plane antiferromagnetic correlations remain spin-liquid at low temperature. Such 
a novel partial spin freezing out of a two-dimensional spin-liquid observed in this cold neutron 
scattering study is likely due to a delicate balance between disorder and quantum fluctuations in 
the quasi- two dimensional 5=1/2 Heisenberg system. 



I. INTRODUCTION 



The parent compound for high transition-temperature 
superconductors, La2Cu04, is an antiferromagnetic in- 
sulator. Magnetic exchange interaction J between the 
nearest neighbor 5=1/2 spins of Cu^"*" ions in the Cu02 
plane is several orders of magnitude stronger than the in- 
terplane exchange interaction, making quantum spin fluc- 
tuations an essential ingredient for magnetic properties in 
the quasi-two-dimensional (2D) Heisenberg systemiiSi^. 
The Neel temperature Tjv of La2Cu04 is suppressed 
rapidly to zero by = 2—3% hole dopants such as 
Sr, Ba or Li^i^, while it is suppressed with isovalent 
Zn substitution at a much higher concentration close to 
the site dilution percolating threshold of ^30%^. The 
strong effect of holes has been shown to be related to in- 
duced magnetic vortices, which are topological defects in 
2D systems^. The paramagnetic phase exposed by hole 
doping at T <IC J /ks is dominated by zero-point quan- 
tum spin fluctuations and is referred to as a quantum 
spin liquid^ . Detailed predictions for spin dynamics have 
been made for the quantum spin-liquid^ ^ 

However, in a wide doping range of La2Cui_a;Li2;04 
below ~10 K, a spin-glass transition has been reported 
in muon spin rotation (/xSR)^, nuclear quadrupole res- 
onance (NQR)^" and magnetizationii studies. A simi- 
lar magnetic phase diagram has also been reported for 
La2-:.Sr^Cu04 and Yi_^Ca:,Ba2Cu30ai2*iiiliSiiLii. 
In conventional spin glasses, magnetic interactions are 
more or less isotropic in space, and the entire spin sys- 
tem is believed to be frozen in the spin-glass phasei^. 
Such was also the conclusion of a comprehensive mag- 
netization study on La2-a;Sra;Cu04^^. Although magne- 
tization can only account for a tiny fraction of spins, 
theoretical pictures were proposed for spin-freezing in 
the whole sampleiLiS. If the spin-glass phase in hole- 
doped cuprates behaved as in conventional spin-glasses, 
the ground state would be a spin-glass, instead of the 
Neel order for doping smaller than Xc, or a quantum spin 



liquid for doping larger than Xc- Thus, as pointed out 
by Hassclmann et al.^", the quantum critical point of the 
antiferromagnetic phase at ~ 2—3% would be pre- 
empted. 

In widely circulating "generic" phase-diagram for lam- 
inar cuprates, the "reentrant" spin-glass transition below 
the Neel temperature is generally ignored. Also generally 
ignored is the spin-glass transition below the supercon- 
ducting transition. The spin glass phase exists side by 
side with the Neel order at lower doping and the super- 
conducting order at higher doping in this neat picture. 
This "generic" picture does not conform to experimental 
results, and serves to support the theory that the spin- 
freezing is an extrinsic dirt effect. However, there are 
other theories which consider spin-freezing intrinsic to 
the doped cupratea^iiSS. Physical quantities in the dop- 
ing regime, including spin excitation spectra, have also 
been calculated from microscopic model^SiSi. 

Recently, 2D spin fluctuations in La2Cui_a:Li2;04 
(0.04 < X < 0.1) were observed to remain liquid-like be- 
low the spin-glass transition temperature^''''^^, ~ 9 K, 
which can be reliably detected using the /iSR technique^. 
The characteristic energy of 2D spin fluctuations satu- 
rates at a finite value below ~50 K^^-— , as expected for 
a quantum spin liquid^, instead of becoming zero at Tg 
as for spin-glass material o^^i^^ . To reconcile these appar- 
ently contradicting experimental results, we have con- 
ducted a thorough magnetic neutron scattering investi- 
gation of La2Cuo.94Lio.06O4 to search for spin-glass be- 
havior. We found that in addition to the liquid-like 2D 
dynamic spin correlations, the rest of spins which par- 
ticipate in almost 3D and quasi-3D correlations become 
frozen in the spin-glass transition. This partial spin freez- 
ing in the laminar cuprate is distinctly different from to- 
tal spin freezing in conventional 3D spin- glass materials. 
The observed phase separation into spin glass and spin 
liquid components of different dimensionality sheds light 
on a long-standing confusion surrounding the magnetic 
ground state in hole-doped cuprates. 

The remaining of the paper is organized as the fol- 
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lows. Section II covers experimental details concerning 
the sample and neutron scattering instrumentation. Sec- 
tion III covers small angle neutron scattering, which is 
the ideal tool to detect ferromagnetic spin clusters pro- 
posed in some theories for the spin-freezing state. Section 
IV covers cold neutron triple-axis measurements. The 
excellent energy resolution is important for this study. 
Finally, in Section V, we discuss and summarize our re- 
sults. 



II. EXPERIMENTAL DETAILS 

The same single crystal sample of La2Cuo.94Lio.06O4 
used in the previous higher energy study^^ was inves- 
tigated in this work. T^siS K was determined in /xSR 
study^ and is consistent with magnetization work^^. The 
lattice parameters of the orthorhombic Cmca unit cell are 
a=5.332l, 6=13.121 and c=5.402l at 15 K. 

Wave-vector transfers q near (000) and (100) in both 
the (hOl) and {hkO) reciprocal planes were investigated 
at NIST using the 30 meter high resolution small angle 
neutron scattering (SANS) instrument at NG7, and cold 
neutron triple-axis spectrometer SPINS. We set the array 
detector of NG7-SANS to 1 and 9 m, corresponding to 
a q range from 0.012 to 0.39 and from 0.0033 to 
0.050 respectively. At SPINS, the (002) reflection of 
pyrolytic graphite was used for both the monochromator 
and analyzer. Horizontal SoUer slits of 80' were placed 
before and after the sample. A cold BeO or Be filter was 
put before the analyzer to eliminate higher order neutron 
in the fixed Ef=S.7 or 5 meV configuration, respectively. 

Sample temperature was controlled by a pumped ^He 
cryostat which could reach down to 1.5 K. 
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FIG. 1: (color) Measured SANS cross-section in the 10 pixels 
wide rectangle shown in the inset at 3 K (squares) and 30 K 
(diamonds), and the difference of the intensity between 3 K 
and 30 K (circles) as a function of wave vector transfer q. 
Lines are guide to the eye. Inset: the SANS pattern at 30 
and 3 K, with the intensity color scale at the left. 



The difference intensity (circles) fiuctuates around zero, 
and its standard deviation sets a upper limit of 1.5x 10^^ 
bn or 1.4xlO~^/iB per Cu for ferromagnetic moments in 
the clusters. 

This result provides serious constrain on the class of 
theoretical models for the spin-glass transition in doped 
cuprates^^ which lead to formation of ferromagnetic clus- 
ters. Instead of this "large spin fixed point", models lead- 
ing to other fixed points such as "Griffiths fixed point" 
as discussed by Lin et alm^ may be considered. 



III. SMALL ANGLE NEUTRON SCATTERING 

Hole induced ferromagnetic exchange has been theoret- 
ically proposed in the Cu02 planei2iS£. It is regarded as 
competing with the original antiferromagnetic exchange, 
thus, leading to the spin-glass transition. Although long- 
range ferromagnetic order has never been observed, there 
is the possibility of short-range ferromagnetic spin clus- 
ters which freeze in the spin-glass state in this class of 
spin-glass models^^'^*. SANS has been demonstrated as 
a powerful tool to probe such clusters^^. 

Two reciprocal zones of La2Cuo.94Lio.06O4 were stud- 
ied, with incident beam parallel to the (001) or (010) 
direction. Therefore, any spin orientation in the sample 
can be detected in our experiment. The experiments were 
carried out at 3, 10, 15, 30 and 80 K. A collection time 
of 1 or 2 hours per temperature provides good statistics. 

No temperature dependence in the scattering patterns 
could be detected. The inset to Fig. ^ shows SANS pat- 
terns at 3 K and 30 K with incident beam parallel to the 
(001) direction. Intensity at 3 and 30 K in the rectangu- 
lar box on the SANS pattern is shown in the main frame. 



IV. TRIPLE-AXIS NEUTRON SCATTERING 

While no appreciable ferromagnetic signal was de- 
tected for La2Cuo.94Lio.06O4, as in other Li-doped 
La2Cu0 4"§f™' , antiferromagnetic scattering was readily 
observed along the rods perpendicular to the Cu02 plane 
and intercepting the plane at the commensurate (tt, tt)- 
type Bragg points of the square lattice. This means 
that antiferromagnetic correlations in the Cu02 plane 
are chessboard-like, which is similar to electron-doped 
La2Cu04^'^^, but different from the more complex, in- 
commensurate ones in La2-2;Srj;Cu04 at similar hole 
doping^^. 

Scans through such a rod in the Cu02 plane at var- 
ious temperatures with the SPINS spectrometer set at 
E—Q are shown in Fig. |2Ia). Inelastic scans have been 
reported previously in a related but different study which 
focuses on scaling in different quantum regimes^^. There 
is little change in the peak width in these scans, consis- 
tent with previous results of temperature independent 
in-plane correlation length for La2Cuo.95Lio.o504i and 
La2_:rSr:rCu04 (0.02< x <0.04)33 below 300 K. Mod- 



3 




0.9 1 
(/lOO) 



(1A;0) 



FIG. 2: Representative magnetic quasielastic and elastic scat- 
tering along a) an in-plane direction and b) the interlayer di- 
rection between 1.5 and 180 K. Open squares in b) were mea- 
sured at (1.06,fc,0) and represent background. The solid lines 
are resolution convoluted {q, E) + S"'^-°(q, E) in Eq. 



eling the width of the rod in Fig. [2Ja) with Lorentzian 



7r[l + {qif 



(1) 



the lower limits from deconvolution is £^2 ^ 274A, where 
the □ indicates the correlation length as in-plane. These 
large antiferromagnetic clusters in the Cu02 plane cor- 
relate in three different ways in the interlayer direction, 
giving rise to almost 3D, quasi-3D and 2D magnetic cor- 
relations. Let us now examine the three components. 

Scans along the rod in the interlayer direction, with 
the SPINS spectrometer set at E—0, were measured 
at various temperatures from 1.5 to 180 K. A few of 
them, at 1.5, 49 and 180 K, respectively, are shown 
in Fig. Elb). Magnetic intensity is composed of both 
broad and sharp peaks at magnetic Bragg points (100) 
and (120) of the parent compound. The (110) peak is 
temperature-independent thus nonmagnetic. Fitting the 
broad peaks to Eq. we obtained an interlayer corre- 
lation length ^'JSf _ 5.2(4) A. Again, no temperature 
dependence can be detected for ^'-^-D }-,giQ-^ 49 }^ Above 
49 K, signal is too weak to have a reliable determination 
of i,'^'^^ ■ Thus, the quasi-3D spin correlations are typi- 
cally three planes thick. For the sharp peak at (100) or 
(120), only the lower limit for the correlation length can 
be reliably estimated: > 168A, since the width is 

close to instrumental resolution. Therefore, the number 
of correlated antiferromagnetic planes is more than 50, 
resembling a 3D antiferromagnetic order. 

Both the broad and sharp peaks in Fig. |2Ib) are 
energy-resolution-limited with the half-width-at-half- 
maximum = 0.07 meV. The energy scan in Fig. |3|is an 
example and more can be found in Fig. 2 in reference 
f2^]. However, these peaks should not be regarded au- 
tomatically as from static magnetic order. Static mag- 



netic signal was observed only below Tg=8 K at the spin 
glass transition in /iSR study^, which has a much bet- 
ter energy resolution. Thus, the quasi-3D and almost 3D 
correlations are slowly dynamic for T > 8 K, with their 
spectra faster than 1 MHz^'^'^, the zero-field /iSR static 
cutoff frequency, but slower than 17 GHz=0.07 meV/h, 
the frequency resolution at spectrometer SPINS. 

The 2D antiferromagnetic correlations have been in- 
vestigated in details^. The dynamic magnetic structure 
factor. 



52^(q,i?) = 5]£fc(/.n)^ 

^ — ' TT 1 



X"{E) 



(2) 



where r is a magnetic Bragg wave- vector and k = q — r, 
has been determined from measurements in the energy 
range, E < 4.2 meV, between 1.5 and 150 K. Eq. Q is 
independent of fc, befitting to a 2D magnetic correlation, 
see the flat k scan at 1.2 meV in Fig. O The almost 
3D and quasi-3D spin correlations described in previous 
paragraphs can be written as 

S'" (q, E) = I'^ J2 ^^"^ ('^n ik-r,)C'/'iE) (3) 

T 

and 

S"^^ (q, E) = I^'^ ('-□)^''" (fc - (E), 

T 

(4) 

respectively, where e < 0.07 meV, the spectrome- 
ter energy resolution. Note that we use conventional 
Lorentzian function, Eq. J^t, to model sharp peaks which 
we could not experimentally resolve, in addition to C^''^ 
in Eq. 10} which we could resolve. We are fully aware that 
the true peak profile can be different for these unresolved 
peaks. The use of Eq. |^ is for the purpose of calculating 
resolution convolution of Eq. l|3l-lEJ, which is used in the 
following paragraphs to obtain correct normalization of 
j2D ^ jqZD j3D r^Yie choice of the function will not 
affect the result as long as the function describes a sharp 
peak significantly narrower than instrument resolution. 

With negligible ferromagnetic correlations, the total 
dynamic structure factor is a summation of Eq. (01-10}, 

S{q, E) = (q, E) + S'^''' (q, i?) + S'^ (q, E). (5) 

Of the four variables of S{q,E), qn are fixed at 
the (tt, 7r)-type Bragg points by the sharply peaked 
(k|-|)34 comprehend the composition of ^(q, E), 
it is sufficient to plot ^(q, E) as a function of E and 
the interlayer wavenumber k. Such a plot of measured 
S{q,E) at 1.5 K is shown with a logarithmic intensity 
scale in Fig. |31 The temperature and q independent in- 
coherent scattering and other background at i? = has 
been subtracted, which can be determined, e.g., by the 
180 K scan in Fig. 0b). The sharp peak fitted by the red 
curve is from S^^ (q, E) , the narrow blue ridge at E—0 
from S"^^^ {q, E), and the green surface from S^^{q,E). 
The red peak at (100) is about one order of magnitude 
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FIG. 3: (color) Measured ^(q, E) as a function of E and in- 
terlayer k at 1.5 K with a logarithmic intensity scale, showing 
three color-coded magnetic components in Eq. (5). S^^{q, E) 
(red) and S"'^°(q,_E;) (blue) are energy-resolution limited at 
-E = 0, representing very slow spin dynamics which is as- 
sociated with the spin-glass freezing. They are modulated 
along the interlayer k direction. The spin-liquid component, 
S^'^{(\,E) (green), has a finite energy scale of about 1 meV 
below 50 K, and O.lSfcBT above 50 It is flat along the 

k direction. A few representative scans at 1.5 K are shown 
with yellow symbols. The black surface indicates background 
of ~ 1.3 counts/min. 



stronger than the peak intensity of the blue surface, and 
three orders of magnitude stronger than the peak inten- 
sity of the green surface. Thus, S'^^{(\,E) is the easiest 
component to be observed in a neutron scattering ex- 
periment, and is often mistakenly attributed to a static 
magnetic order. 

The spectral weights J dqdE S^^ {q, E)=I^° and 
f dqdE S''^^ {q, E)=I'}^^ can be obtained by fitting 
resolution-convoluted Eq. Q-lB to scans such as those 
shown in Fig.|2Ib). They are shown as a function of tem- 
perature in Fig. ^ with I^^ magnified by a factor of 5 
for clarity. For the 2D component, the spectral weight is 



= / dE 



2x"iE) 



(6) 



where the integration limits are ±cxi. Green squares in 
Fig. 21 represent the lower bound of with the en- 
ergy integration limited in \E\ < 10 meV, using the 
analytical expression of x"{E) in reference f2^] to ex- 
trapolate to E—IO meV, where spin fluctuations were 
observed in La2Cuo.9Lio.1O4 using a thermal neutron 
spectrometer—. 

and /'■^^ appear simultaneously below ^150 K. 
Their concave shape in Fig. 0] differ drastically from the 
usual convex-shape of a squared order parameter, orange 
circles, which was observed in ^SR study below Tg=8 
They are typical neutron scattering signal from slow dy- 
namic spin correlations in spin-glasses^^i^, which fluc- 
tuate in the frequency window between 1 MHz and 17 



FIG. 4: (color) Temperature dependence of spectral weights 
I^^ (green), /'■^^ (blue) and I^^ (red) in the same unit for 
three experimentally separable antiferromagnetic components 
in Eq. @ . + is the total spectral weight of the spin- 
glass component. I^^ is the spectral weight of the spin-liquid 
component within l^l < 10 meV, thus the lower limit of its 
total spectral weight. The orange circles represent squared 
static order parameter of the spin-glass transition, which was 
measured by ^iSR^ and equals to + 1'^^° at T=0. 



GHz for T > 8 K, and below 1 MHz for T < 8 K. Previ- 
ously, energy-resolution-limited neutron scattering from 
Lai.94Sro.o6Cu04 was observed to have a similar temper- 
ature dependence as I^^ in Fig. ^ and was attributed 
to spin freezingi^. The kink of at 20 K reflects an 
increased Tg from 8 K to 20 K when probing frequency is 



increased from 1 MHz to 17 GH: 



,13,18 



At Hz, r„ « 6 K 



from DC magnetization measurementsii. The increase of 
Tg with measurement frequency is a hallmark of glassy 
systemai^. 



V. DISCUSSIONS AND SUMMARY 

The fact that /^^ decreases below rg(17GHz)«20 K 
while I'^^^ continues to increase indicates that the 
"Edwards- Anderson order parameter"i2i2L2^ distributes 
only along lines such as the (l/cO). In conventional 
spin-glasses, the "Edwards-Anderson order parameter" 
is more isotropically distributed in the q-spaceiSiSLS. 
Thus, the spin-glass state in La2Cuo.94Lio.06O4 is charac- 
terized mainly by interlayer disorder which upsets phase 
correlation between large antiferromagnetic clusters in 
different Cu02 planes. This picture offers a possible 
alternative to the conventional competing antiferromag- 
netic/ferromagnetic interaction model for spin freezing in 
doped cuprates. In addition, it suggests that the weak 
interlayer exchange interaction likely plays an important 
role in the finite temperature spin-glass transition in the 
quasi-2D Heisenberg magnetic systems. 

Another important difference from conventional spin- 
glasses in which all spins are believe to freeze at low 
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temperature is that only a fraction of spins freeze in 
La2Cuo.94Lio.06O4. Other spins in 2D correlations re- 
main fluctuating down to 1.5 K. This is consistent with 
numerical evidence that quantum fluctuations prevent 
spin-glass transition for 2D 5*= 1/2 Heiscnberg system^S;. 
The spin-glass component in our sample has to acquire 
interlayer correlations to achieve a higher dimension in 
order to be realized. It appears that the lower critical 
dimension for a 5=1/2 Heisenberg quantum spin glass is 
between 2 and 3. 

A further difference from conventional spin-glasses, 
for which one can measure the narrowing of magnetic 
spectrum toward E=0^, is that when S^^ (q, E) and 
S''^^{c[,E) in La2Cuo.94Lio.06O4 become detectable at 
about 150 K, they are already energy-resolution-limited, 
with spins fluctuating much slower than 17 GHz. This 
property of S^^{q,E) and S"?^^(q, i?) resembles the 
classic central peak phenomenon in the soft phonon 
transition'^^-'^* . The disparate dynamics of the central 
peak and phonon are explained by Halperin and Varma™ 
using a phase separation model: defect cells contribute 
to the slow relaxing central peak while coherent lattice 
motions (phonons) to the resolved inelastic channel. This 
mechanism has been applied with success to a wide class 
of disordered relaxor ferroelectrics^Si^. 

For La2Cuo.94Lio.06O4, we envision that disorder ac- 
companying doping prevents the long-range order of 
the antiferromagnetic phase mainly by upsetting inter- 
layer magnetic phase coherence, see Fig. O for the q- 
distribution of frozen spins. This upsetting is not uni- 
form in the Grifflths fashion^ with weak and strong cou- 
pling parts in the sample. In our laminar material, how- 
ever, the weak and strong coupling parts have different 
dimensionality: 2D and nearly 3D, respectively. The 2D 
part is a spin liquid and represents essentially the whole 
system at high temperature, see Fig. 0] Part of sample 
with stronger interplane coupling tends to order three 
dimensionally below ~150 K, producing S^^{q,E) and 
593^ (q,£;). The condensat ion of the 2D spin liquid at 
^ 150 K into the quasi-3D dynamic clusters of dimin- 
ishing energy scale, instead of a true long-range order, 
may reflect the divergent fluctuations which destabilize 
static order at finite temperature for 2D random XY or 
Heisenberg systema^SiM^i. The nearly 3D spin-glass in- 
stead of a 3D antiferromagnet finally orders at a much 



reduced Tg « 20 K, when /^^^ + approaches the 
2D spectral weight (Fig.^. The coexistence of spin hq- 
uid and spin glass components at low temperatures may 
be a general consequence of no "mobility edge" separat- 
ing finite and infinite range correlations for a 2D random 
system^. Recently, Monte-Carlo simulations of a doped 
2D classical antiferromagnet suggest that there are two 
populations of spins: one with fast and the other with 
slow dynamics^. This is consistent with our experimen- 
tal results and the Griffiths picture for random magnetic 
systems. A phenomenological Halperin and Varma model 
may be built for spin dynamics in doped cuprates based 
on these microscopic insights. 

In summary, spins in La2Cuo.94Lio.06O4 develop dy- 
namic antiferromagnetic order in the Cu02 plane with 
very long below 180 K. The characteristic energy of 
the 2D spin fluctuations is O.lSfcsT for T > 50 K and 
1 meV for T < 50 }<M. Below --150 K, interlayer phase 
coherence appears between some of these planar antifer- 
romagnetic clusters with an energy scale smaller than 
70 /ieV. While the 2D antiferromagnetic correlations in 
an individual plane remain liquid down to 1.5 K, co- 
herent multiplane antiferromagnetic correlations become 
frozen below Tg. The phase separation into 2D spin- 
liquid and spin-glass of higher dimension with an unusual 
q-structure for the "Edwards- Anderson order parameter" 
is most likely related to quasi-2D nature of magnetic ex- 
change in the cuprates and is distinctly different from 
conventional spin-glasses. 

A theory of spin-glass in doped cuprates should include 
interlayer coupling. Theory explaining both the partial 
spin freezing and the observed crossoveriS^sS^ of quantum 
spin fluctuations are called for. The heterogeneous mag- 
netic correlations, instead of a uniform magnetic phase, 
suggests the possibility that superconductivity and the 
almost 3D antiferromagnetic order may reside in differ- 
ent phases in La2-a;Sra;Cu04 and Yi_a;Caa;Ba2Cu306+j,. 
Similar, detailed q, E and T dependent cold neutron 
spectroscopic study on these cuprates are desirable. 
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